Transverse radiation force in a tailored optical fiber 
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We show, by means of simple model calculations, how a weak laser beam sent though an optical 
fiber exerts a transverse radiation force if there is an azimuthal asymmetry present in the fiber such 
that one side has a slightly different refractive index than the other. The refractive index An needs 
only to be of very small, of order 10 -3 , in order to produce an appreciable transverse displacement 
of order 10 fim.. We argue that the effect has probably already been seen in a recent experiment 
of She et al. [Phys. Rev. Lett. 101, 243601 (2008)], and we discuss correspondence between these 
observations and the theory presented. The effect could be used to bend optical fibers in a predictable 
and controlled manner and we propose that it could be useful for micron-scale devices. 



I. INTRODUCTION 

Recent years have seen an increased interest in radia- 
tion forces in optics. Optical tweezers, atom traps, opti- 
cal manipulation of soft materials such as interfaces be- 
tween liquids - especially near the critical point where 
surface tension is small - are typical examples. This trend 
will in all probability continue in the time to come. 

Our objective in the present Communication is to point 
towards the possibility - to our knowledge not so far con- 
templated in the literature - to create a tailored trans- 
verse optical force in a fiber transmitted by a laser beam. 
The beam may be pulsed, or it may be continuous. The 
clue of the principle is to introduce an accurate mechan- 
ical imbalance in the fiber, implying a slight asymmetry 
in the refractive index n. (In practice, such a deviation 
from axisymmetry may easily result inadvertently, dur- 
ing the mechanical drawing of the fiber.) If one side of 
the fiber is harder than the other, there may be a slight 
refractive index difference An between the two sides, re- 
sulting in a transverse optical force. As fibers of micron 
scale cross sections are very light and bend easily, a side- 
ways motion may easily occur. The effect, besides being 
of basic interest, may be of practical utility. We will de- 
scribe the effect, making use of simple models for the 
fiber, and thereafter compare with a recent experiment 
which in our opinion has most likely already observed 
this effect. 

The problem is to some extent related to the one- 
hundred years old Abraham- Minkowski debate on the 
correct electromagnetic energy-momentum tensor in di- 
electric media. From a physical point of view the key 
issue is that one is dealing with a nonclosed system, mat- 
ter and field. Macroscopic of phenomenological electro- 
magnetic theory implying the use of a permittivity e and 
permeability \x means that one is dealing with a com- 
plicated interaction system involving external fields, in- 
ternal fields, and constituent molecules, by using only 
simple material parameters. The solution to the prob- 
lem lies in extracting the energy-momentum form that 
leads to a theoretical description of observable effects 
in a clean and simple manner. For an overview of 
the Abraham-Minkowski debate, cf. e.g. [l[. There ex- 



ist at present a great number of papers discussing the 
Abraham-Minkowski problem, for instance the recent 
Ref. and the review 

Assume now for definitcness that the fiber is vertically 
hanging, and that a laser pulse is transmitted through 
it. The general expression for the electromagnetic force 
density in the medium, assumed hereafter nonmagnetic, 
is derived from a given stress tensor as /j = dk^ik- 
For the Abraham and Minkowski tensors this implies (cf. 
Refs. d Q for details) 
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(1.1) 



Here f = — (sq/2)E Vu is nonvanishing in any re- 
gion where n varies, especially in the surface regions. 
We use notation E(r,t) = Re{E(x)e i(wt - /32) } and skip 
the exponential factor in the following (cf. Q Ch. 1.3, 
for notational details). As this force is common for the 
Abraham and Minkowski tensors, it may appropriately 
be denoted as f AM . The Abraham momentum density 
g A = (l/c 2 )E x H occurs in the second term. We will 
henceforth ignore it, due to the following two reasons: (i) 
for a stationary beam the term fluctuates out when aver- 
aged over an optical period; (ii) under perfectly axisym- 
metric conditions the force exerted on the fiber during 
the transient entrance and exit periods has necessarily 
to be vertical, thus being unable to initiate any sideways 
motion. 

In the following we investigate the effect of the force 
f AM when the refractive index contrast An is assumed 
known from the mechanical production process. Two 
simple planar models will be considered, in order of in- 
creasing complexity. Finally we compare the theory with 
a recent experiment by She et al. [6| and demonstrate 
how the deflection observed could very plausibly be a 
demonstration of the effect mentioned. 



II. SLAB OPTICAL WAVEGUIDE 

We will henceforth limit ourselves to planar geometries 
for the fiber. This is mathematically simplifying, but the 
model is nevertheless expected to incorporate the essen- 
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tials of the imbalance effect. Probably the simplest ar- 
rangement is to consider a uniform slab, infinite in the 
horizontal y direction, having a finite width 2a in the 
other horizontal x direction. 
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FIG. 1: Left: The simplest model, a slab waveguide with a 
slight difference in refractive index above and below. Right: 
The geometry of a vertically hanging optical fiber subject to 
sideways motion. 



The setup is sketched in Fig. [Th; it is essentially the 
same as Fig. 2.1 in Ref. The beam is propagating 
into the plane, in the z direction. On the lower side 
of the slab [x < —a) we assume vacuum (or air), with 
refractive index no = 1. On the upper side {x > a), we 
assume that there is a dilute medium, extending to x — 
oo, with refractive index ri2 = 1 + An, An <C 1. When 
light propagates through the fiber, there will thus be an 
imbalance in the surface force densities on the lower and 
upper surfaces. Let use choose a TE mode (following the 
notation in Ref. @). It means that the electric field has 
only an y component different from zero, E = (0, E y ,0). 
fn the dielectric boundary layers located around x = —a 
and x = a, the transverse force component is f x M — 
— (£q/2)E 2 d(n 2 )/dx, yielding the respective surface force 
densities to be 

a x (-a)=-le Q E 2 y (-a)(n 2 -l), (2.1) 
a x (a)=-\e Q E 2 y {a){nl-n 2 ). (2.2) 

We have here taken into account that the longitudinal 
component E y is continuous across the surfaces, and that 
(E 2 ) = \E 2 y where ( ) is average over oscillations in time 
and the z direction. 

The net transverse surface force density is a x = 
+ Ox(tt)- As An is small, we can make use of 
the expression for E 2 corresponding to a symmetric fiber, 
thus with the assumption ri2 = 1. We get accordingly to 
first order in An 



cr x = -±e Q E 2 (a) ■ An, 



(2.3) 



where now E 2 (a) refers to the symmetric situation. We 
see that for positive An, the surface force is directed 
downwards in Fig. [TJi, i.e., in the negative x direction. 
This is as it should be, as surface forces are always di- 
rected towards the optically thinner region at a dielectric 
surface. 



The presence of a x makes it possible to regard the fiber 
as an elastic rod exposed to a constant transverse load. 
Let us first, however, relate a x to the total power P in 
the fiber. For the TE mode we have, using the same 
notation as in Ref. [5|, 



A cos(na 
E y = I A cos(kx — 0), 



,4cos(Ka + (j))e^ x+a \ 



x > a 

—a < x < a 
x < —a. 



(2.4) 



Here k = v/n 2 u? 2 /c 2 — /3 2 and £ = \//3 2 — ui 2 /c 2 , where 
/3 lying in the interval w/c < (3 < ntojc is the wave 
number component in the z direction. The corresponding 
nondimensional transverse wave vectors are u = kci, w = 

The electromagnetic boundary conditions, requiring 
that dEy/dx be continuous, yield the following equations: 



u = ^rmr + arctan (w/u) , cf> = ^rmr, 
with m = 0, 1, 2..., while u and w are related via 



to 2 a 2 



w 
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(2.5) 



(2.6) 



From the above equations (3 and w can be calculated, and 
we can find the relationship between P and the constant 
A in Eq. (O using formula (2.34) in Ref. 0|: 



A = 



2w/ioP 



0ab{l + l/w) 



(2.7) 



Recall that P refers to the total power transmitted by 
the fiber. In the planar model, we have let b denote the 
fiber width in the y direction. The cross sectional area 
of the model fiber is thus 2a6, and the power per unit 
length in the y direction is P/b. Edge effects because of 
the finite value of b are ignored. The transverse surface 
force density can now be expressed as 



abn e c(l + l/w) n 2 — 1 



■At 



where we have defined 



Pc/w 



(2.8) 



(2.9) 



and used that the continuity of dE y /dx across the inter- 
face at x — a implies that 



cos 2 (m- 



cos 2 (u + < 



i 



(2.10) 



For practical purposes it may be convenient to express 
An in Eq. (|2.8p in terms of the corresponding increase 
Ap in material density. We can make use of the Clausius- 
Mossotti relation, which is a good approximation at least 
for nonpolar materials. We then get 



An 2 



1 

3P 



(n 2 - l)(n 2 + 2) ■ Ap. 



(2.11) 
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Consider next the fiber as an elastic rod of rectangular 
cross section, clamped at one end (z — 0) and free at the 
other end (z = L). For convenience we let the z axis be 
horizontal. We choose b — 2a, implying a square cross 
section of the fiber. The transverse load per unit length in 
the longitudinal direction is a x b acting downwards in the 
negative x direction. The governing equation for the elas- 
tic deflection is (we ignore gravity) x""(z) — a x bj(El), 
where E is Young's modulus and / = 6 4 /12 the moment 
of inertia of the cross-sectional area about its centroidal 
axis Q ■ The solution of the governing equation is 



x{z) = - 



<J x b 



l (z 2 — ALz + 6L 2 



(2.12) 



The deflection at the tip, called simply Ax, is thus 



Ax 



i x bL 
8EI 



(2.13) 



We can now calculate the perturbation An required to 
produce a relative deflection Ax/ L at the tip: 



lAnl 



8EI 
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n„ac 1 
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\Ax\ 



(2.14) 



In order to obtain an order of magnitude estimate for 
the required difference in n in order to yield a prescribed 
value for Ax we insert numerical values that are appropri- 
ate for a low-intensity laser beam in a fiber (the numbers 
are typical, and the same as used in Ref. Q): A = 650nm, 
P = 6.4mW, L = 1.5mm, b = 2a = 450nm. 

For the eigenvalues of the transverse wave number 
(3 (or, equivalently, n e ), corresponding to the guiding 
modes of the planar fiber, we use again the solutions 
for the symmetrical situations, since corrections to this 
only enter beyond leading order in An. Following [f|, the 
eigenvalue of index m = 0, 1, 2, ... solves 



fib) = rv J - b - arctan \/ -^—^ - ™ = 0, (2.15) 



wherew = *fVn 2 - 1 and b = (n 2 -l)/(n 2 -l). With the 
current numbers there are only two modes, n ej o = 1.3635 
and n e4 = 1.1072. 

As an example, assume now we desire a lateral dis- 
placement Ax = 15 urn. With the above numbers we 
find that the required An is only 5.5 • 10~ 4 for m = 
and 1.6 • 10~ 4 for m = 1. Minor mechanical defects from 
production could easily give rise to changes in the reflec- 
tive index of this magnitude. While the planar model 
is likely to underestimate the required An slightly com- 
pared to a circular fiber, as an order of magnitude esti- 
mate it demonstrates the feasibility of the scheme. 



III. FOUR-LAYERED MODEL 

Regard now a slightly more realistic model where the 
slab has a layer of slightly higher refractive index on one 



(3) Vacuum 
(2) — n+An 
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FIG. 2: A two-layered model. 



side. The geometry is as considered in Fig. [2] where the 
slab, still of width 2a and refractive index n, has refrac- 
tive index raised by An in a layer of thickness Aa. Again 
we consider the TE mode, whose solution for the electric 
field component E y is written [in layers (3) to (0) from 
top to bottom] 



E„ = A 



cos(u — Au — 4>) cos(Aw — tp)e~^ x ~ a ^ 
cos(u — Au — (f>) cos(k'(x — a + Aa) — ip), 
cos(kx — 4>) cosi/' 
cos(u + 4>) cos i/j ■ et( x+a ) 



where Au = kAa, At; 



'Aa, 



(3.1) 

and ?p are unknown 



phase angles and n' 2 = (n + An) 2 k 2 — (3 2 « K 2 + 2nfc 2 An. 
The net surface force per unit area is now given by a x = 
E 2 (x)\:Z a . 

a x = - \e a A 2 [cos 2 (u - Au - 4>) cos 2 (Aw - -0) 

-cos 2 (u + </>)cos 2 -0]- (3.2) 

To leading order in An we may again use the relation 
(12. 7|) for A, and cos 2 (u + 4>) is again given by equation 
(l2~T0l) . 

From continuity of dE y /dx at x = a it follows that 
cos 2 (Aw - ijj) = k' 2 /(£ 2 + k' 2 ), that is, 



cos 2 (Aw — ijj) sa cos 2 (u + < 



1 - 



2nAn 
n 2 - 1 



2nAn 
n 2 - 1' 



Similarly, from the condition of continuity of dE y /dx at 
x = a — Aa we derive that, to linear order in An we have 



cos 2 (u — Au 



cos V 1 



2nAn 



• sin ip 



Thus we find with some calculation that the surface force 
per unit area may be written as 



heoA 2 



lAr 



1 



cos ip 



1 



cos 2 ?/> 



(3.3) 



We may express cos 2 ip by means of the parameter Au 
by using the equation of continuity of dE y /dx at x = a, 
which may be written ip = — arctan(£/K) + Aw — mir, 
m = 0, 1, to derive the relation, valid to leading order 
in An: 



cos 2 ip 



n e ( 

tI 



cos Au - 



sin Au 



(3.4) 
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FIG. 3: (a) Required An as a function of Aa/a to obtain 
Ax = 15[J.m with parameters as given in text. The thin 
straight lines show the corresponding values for the simple 
model of figure [T] (b) Surface force per unit area divided by 
An as function of Aa/a. 



When An is small we may use the same eigenvalues of 
n e as before. 

Assuming once again that a displacement of Ax — 
15u.m is sought, figure [3^, shows how the required An 
changes with varying values of Aa/a. For the m = 
mode somewhat higher values of An when Aa/a is small. 
This is because a x oc AnAa as Aa — > 0. The required An 
remains small, however, as long as the layer of slightly 
increased n is non-vanishing. For An > 1% for example, 
Aa/a must be smaller than about 0.035 to achieve the 
observed deflection. 

The surface force density a x / An to leading order is 
plotted as a function of Aa/a in Fig.[3jD. As in the sim- 
plest model of Fig.[T]the force is directed downward when 
Aa/a is smaller than approximately unity. In the present 
example it is reasonable to assume that it is a small num- 
ber compared to 1. We have performed the calculation 
also using different laser wavelengths, in which case the 
graph of Fig. looks qualitatively similar, but the ab- 
scissa is rescaled. Since the point where the force density 
changes sign (about Aa/a = 0.75 in Fig.Oa) then moves 
along the abscissa, it is possible to change the sign of the 
force by changing the optical frequency, provided Aa/a 
is chosen carefully. 



IV. COMPARISON WITH AN EXPERIMENT. 
CONCLUSIONS AND OUTLOOK 

The recent experiment of She et al. [6| appears to 
be a natural example of application for the present the- 
ory. This experiment tested precisely the sideways mo- 
tion of a vertical fiber when transmitted by a laser beam. 
Whereas the authors actually related their observations 
to the Abraham-Minkowski problem mentioned above, 
we do no think that such a conclusion is right, and one of 
us has discussed this issue in more detail elsewhere [8[ . In 
the experiment @ one side of the fiber was slightly harder 
than the other from fabrication Q, and it is natural to 
suggest that this could have given rise to an appreciable 
An. 

Noting how even very small An can cause deflections 
of the magnitude observed, and how the surface force de- 
pends on the geometry of the hardened layer in a nontriv- 
ial way, it is much more likely in our opinion that all the 
observed transverse deviations reported in [|| were due 
to the force f AM . As such the experiment is a striking 
demonstration of the theory presented herein. To ease 
comparison, our numerical examples above used data 
taken from the experiment including a deflection in 
the order of 15|i.m. A change in the refractive index on 
one side of the fiber in the order of less than 1% is shown 
to be sufficient to produce deflections of such magnitude. 

In conclusion, only a minute change in the refractive 
index, around 10~ 3 , in a layer on on one side of an optical 
fiber of centimetre length and micron scale cross-section 
is sufficient to produce an appreciable transverse deflec- 
tion of the fibre. A tailored fiber of this kind, where An 
is prescribed accurately, should provide the possibility 
to bend fibers in a controlled and predictable way. Be- 
cause the lateral force acting on the fiber depends on the 
laser frequency in a non-trivial way, it would moreover 
be possible to design a fiber which bends one way at one 
frequency, and changes direction of movement when a dif- 
ferent frequency is used. The force can be used to carry 
a small load and one can envisage uses within micro- and 
nano robotics. 
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